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Abstract 

> : 

The coherent state of the nonlinear oscillator having nonlinear spectrum, is constructed 
using the Gazeau-Klauder formalism. The weighting distribution, the Mandel parameter are 
studied. The details of the revival structure arising from different time scales underlying 
the quadratic energy spectrum are investigated by the phase analysis of the autocorrelation 
O ' function. 

> : 

• i-H , 

1 Introduction 

The nonlinear differential equation 

(1 + Xx 2 )x - {Xx)x 2 + o?x = 0, A>0 (1) 

was studied by Mathews and Lakshmanan in [1,2] as an example of a non-linear oscillator 
and it was shown that the solution of (1) is 

x = Asin(ujt + 0) (2) 
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with the following additional restriction linking frequency and amplitude 

It has been shown in [1] that Eqn.(l) is obtainable from the Lagrangian density 

' 1 '-2 „.2„2 



- 2(l + A^- aV > < 4 > 
Recently in a series of papers [3-5] this particular nonlinear system has been generalized to 
the higher dimensions and various properties of this system have been studied thoroughly. In 
ref [6], the Schrodinger equatiion corresponding to this nonlinear oscillator has been solved 
exactly as a Sturm-Liouville problem and A-dependent eigenvalues and eigenfunctions were 
obtained for both A > and A < 0. It is to be noted that 

1) this A-dependent system can be considered as a deformation of the standard harmonic 
oscillator in the sense that for A — > all the characteristics of the linear oscillator are 
recovered. 

2) in the Lagrangian (4), the parameter A is present not only in the potential ( i+\ x 2 ) but 
also in the kinetic term. So this nonlinear oscillator may also be considered as a particular 
system with position dependent mass m(x) = ^ 1+ \ x 2^ ■ 

Schrodinger equation with a position dependent mass has found applications in the field 
of material science and condensed matter physics, such as semiconductors [7], quantum wells 
and quantum dots [8,9], 3He clusters [10], graded alloys and semiconductor heterostructures 
[11-19] etc. It has also been found that such equations appear in very different areas. For 
example, it has been shown that constant mass Schrodinger equation in curved space and 
those based on deformed commutation relations can be interpreted in terms of position 
dependent mass [20,21]. This has generated a lot of interest in this field and during the past 
few years various theoretical aspects of position dependent mass Schrodinger equation have 
been studied widely [22-39]. 

On the other hand coherent states have attracted considerable attention in the literature 
[40,41]. Coherent states are generally constructed by (i) using the displacement operator 
technique or defining them as (ii) mininum uncertainty states or (iii) annihilation operator 
eigenstates. However, even when such operators do not exist, different approaches [42-47] 
have been utilized to construct coherent states corresponding to different quantum mechan- 
ical potentials [48-74]. Coherent states of systems possessing nonlinear energy spectra are 
of particular interest as their temporal evolution can lead to revival and fractional revival, 
leading to Schrodinger cat and cat-like states. For constant mass Schrodinger equation, 
potentials like, Poschl- Teller, Morse and Rosen-Morse lead to nonlinear spectra. Time evo- 
lution of the coherent states for these potentials [75-83] is a subject of considerable current 
interest as they can produce Schrodinger cat and cat-like states. 
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In this paper we shall construct coherent states for the nonlinear oscillator (1) and study 
their revival dynamics. The motivation comes from the fact that the study of the temporal 
evolution of a free wave-packet with position dependent mass inside an infinite well [84] has 
revealed that revival and partial revivals are not only different from the constant mass case 
but also they are very much dependent on the mass function m(x). So it will be interesting 
to study the role of the mass parameter A of the nonlinear oscillator (1) on the temporal 
evolution of it's coherent state. 

The paper is organized as follows. In section 2. the Gazeau-Klauder coherent state 
[45,46] for the nonlinear oscillator is constructed and is shown to satisfy the conditions such 
as, continuity of labelling, resolution of unity, temporal stability and action identity. The 
revival dynamics of the coherent state is studied in section 3. In section 4. we summarize 
our results. 



2 Gazeau-Klauder coherent state for nonlinear oscilla- 
tor 

The position dependent mass Schrodinger Hamiltonian H for the nonlinear oscillator(l) is 
given by 



H 



h 2 , , 2 , d 2 h 2 d l / x 2 

— 7^ — (1 + Xx -— — - — Ax-; — h -q — ;r 

2m v ' dx 2 2m dx 2 y \l + A:r 2 



(5) 



The parameter A may be positive or negative. However for A > the discrete energy spec- 
trum can be shown to be finite and consequently for completeness property the continuum 
has to be taken into account. On the other hand for A < there are only discrete energy 
states and henceforth we shall consider this choice. Replacing the parameter g in Eqn.(5) 
by ma 2 + Xha as was done in ref[6], the Schrodinger equation is given by 



-h 2 . 2 , d 2 h 2 d 1 h\ 

■Xx 



2m (1 + XX ] dx 2 ~ 2m 



x 



H — ma(a H )( 

dx 2 v m JK l + Xx 2 



(6) 



Under the transformation x = y -^y ; A = ^ A, the above equation transforms to 

d 2 . d 



+ Ay— -(1 + A)- 



y 



dy 2 a dy v ' (1 + Ay 2 ) 



+ 2e 



V> = 



(7) 



where E = e(ha). For A < 0, the eigenvalues e n and the eigenfunctions ip n of the above 
transformed Hamiltonian are given by [6] 



i> n = (1 - \A\y 2 )Tm)H n (y, A), e n = (n + -) + -n 2 |A|, n = 0, 1, 2 • • -m ■ 



(8) 
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where H n (y,A) are called A-deformed Hermite polynomial [6]. So the eigenvalues of the 

n{n + p) 



Hamiltonian H 1 = H — are 



El = E n -^ = ha 



n + -n 2 |A| 
2 



ha- 



= haei 



// 



(9) 



where p — The Gazeau-Klauder coherent state [45,46] for this system is given by 



|J,7>= 



E 



(J)2ea;p(-?7e^) 



where 7 = at. The normalisation constant N(J) is given by 



n > 



(10) 



iV(J) 



7n 



, < J < R = lim sup n -, +00 pZ 



(11) 



where i? denotes the radius of convergence and p n denotes the moments of a probability 
distribution p(x) 



pn= x n p(x)dx = TTe- ; po = 1 
■'o i=i 



(12) 



For the coherent state (10) 



^ji(i + p) r(n + l)r(n + l + /i) 

11 ~ - ..„-n/t , >P0-1 



i=l 



so that R is infinite and 



/i«r(i + /i) 

2/i(J/i)'' 



r(i + / i) 

K u (cx) being modified Bessel function [85]. 
Also 



(J/i)2 



(13) 



(14) 



(15) 



where I^cx) is the modified Bessel function [85]. 
So the coherent state (10) finally becomes 



v TtTT^) y, (J / ,)fe-^ + ^ 



n > 



(16) 



Below we shall see that the coherent state (10) satisfies the following conditions: 
1. Continuity of labelling: From the definition (10) it is obvious that 
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{J'1)^{J,l)^\J',1 >^|J,7> 
2. Resolution of unity: 

j | J, 7 >< j, 7 |d^(j, 7 ) = j V,7 >< J,iW{J,i) = ^J d ~i J k (J)\J,i >< Jn\dJ 

(17) 

where fc(J) is defined by 



fc(J) = iV(J) 2 p(</) > o, < J < i? 
= p(J) = 0,J>R 



iH L t , 



(18) 



For the coherent state (10) 

k(J) = 2»Ip(2y/jil)K ll (2y/jjl) (19) 

so that the resolution of unity is satisfied 

J \J, 1 xJ, 1 \d f i(J, 1 ) = l (20) 



3. Temporal stability: 

e~^\ J, 7 >= | J, 7 + at > (21) 



4. Action identity: 

< J^ifV.T >= MJ (22) 
The overlapping of two coherent states is 



!> + !) - (jjy)fe-(v-^ 

<J ' 7|J ' 7>= iV(J)iV(J0^ n!r(n+l + ^) (23) 



If 7 = 7 ', the overlapping is reduced to 

^'^^ Tmdmm '^"'^ (24) 

3 Revival dynamics 

In this section we shall study the revival dynamics of the coherent state (10). To demonstrate 
the role of the mass parameter A on the revival dynamics, all the figures below are drawn for 
a fixed value of J = 10 and two different values of /i which is inversely proportional (/i = j|j) 
to the mass parameter A. For a general wave packet of the form 

m)>=J2 c " e ~ iEnt/H \ n> ( 25 ) 

n>0 
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with J2 n >o l c «l 2 = 1' concept of revival arises from the weighting probabilities |c n | 2 . For 
the coherent state (10), the weighting distribution is given by 

( J "»" + * (26) 



nW(n + 1 + fil^y/Jji) 



Since the weighting distribution \c n \ 2 is crucial for understanding the temporal behavior of 
the coherent state (10), we show the curves of \c n \ 2 for J = 10 and different \i in Figure 1. 




5 10 15 20 



Figure 1: Plot of the weighted distribution |c n | 2 given in (26) for J=10 and ji = 28,80. 



The Mandel parameter Q is defined by 



Q = 



(An) 2 

< n > 



(27) 



where 



< n > 
An 



E 

n=0 



lJ n 



N(jy Pri 



(28) 



[< n 2 > - < n > 2 ]5 



The Mandel parameter determines the nature of the weighting distribution function \c n 



The case of Q = coincides with the Poissonian weighting distribution 



< n > n e- <n> 



while 



the cases of Q > and Q < correspond to the super-Poissonian or sub-Poissonian statistics 
respectively. For the coherent state (10), 



Q = v J/i 



where 



< n > 



y 2 (2y/J^) T^igygg 
J M+I (2v/J^) I^yfTJl) 



(29) 



(30) 



6 



In Figure 2 we plot the Mandel parameter Q for J = 10 and [i = 28, 80. It is evident from 
the figure that the Mandel parameter is sub-Poissonian and it has been observed that it 
remains so for all values /i. 




5 10 15 20 25 

J -> 



Figure 2: Plot of the Mandel parameter Q given in (29) for J = 10 and fi = 28, 80. 

Now assuming that the spread An = [< n 2 > — < n > 2 ]^ is small compared to < n >^ n, 
we expand the energy E\ in a Taylor series in n around the centrally excited value n: 

K*>Ek + E l n{n - n) + („ - nf + (n - nf + • • • (31) 

where each prime on E\ denotes a derivative with respect to n. These derivatives define 

2-nh 2nh 
distinct time scales, namely the classical period T d = , , the revival time t rev = — — prr and 

I | 2 I Efi I 

so on. For El in Eqn.(9), T c i = — -— and t rev = — — . There is no superrevival time here 

a(2n + n) a 

because the energy is a quadratic function in n. It is convenient to describe the wavepacket 
dynamics by an autocorrelation function 

A(t) =< 4>(x,0)\i>(x,t) >= |c„|V^ (32) 

n>0 

For the coherent state (10) the autocorrelation function (32) is given by 

A{t) =< J,0\J,at >=^#E ^/^r r e-*** 1 ^ (33) 

N(jy ^n!r(n + l + /x) 

From (33) it follows that \A(t — t rev )\ 2 = \A(t)\ 2 so that the autocorrelation function is 
symmetric about ^ L f L . In other words whatever happens in [0, is repeated subsequently. 
We note that < |A(t)| 2 < 1 and it denotes the overlap between the coherent state at t — 
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and the one at time t. So, the larger the value of |A(t)| 2 , the greater the coherent state 
at time t will resemble the initial one. In Figure 3(a), 3(b) and 3(c) we plot the squared 
modulus of the autocorrelation function (33) in units of t rev for J = 10 and /i = 1,28,80 
respectively. Figure 3(a) clearly shows the full revival. The sharp peaks in Figure 3(b) and 
3(c) arise due to fractional revivals. 
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Figure 3: Plot of the square modulus of the autocorrelation function |A(t)| 2 given in (33) 
for (a) J = 10, n = 1, (b)J = 10, /i = 28, (c) J = 10, /i = 80. 

We can see from figures 3(a) to 3(c) that different values of /x lead to qualitatively different 
types of motion of the coherent wavepacket. 

Now we shall study the mechanism of the fractional revival by phase analysis [86,87]. For 
this, we use the time scale determined by the period of the complete revival t rev = ^ = I. 
So the phase of the n-th stationary state is 

(f) n (t) = 2ir(nn + n 2 )t (34) 

where /i — y|| is assumed to be an integer. 

At arbitrary moments of time, the phases (34) of individual components of the packet are 
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uniformly mixed so that it is not possible to make any definite conclusion about the value 
of the autocorrelation function. However, at specific moments, the distribution of phases 
can gain some order; for example, the phases can split into several groups of nearly constant 
values. The fractional revival of q-th order is defined as the time interval during which the 
phases are distributed among q groups of nearly constant values. 

To consider fractional revival of order q, in the vicinity of time t = ^, it is convenient to write 
n as n = kq + A, where k — 0, 1, 2, • • • and A = 0, 1, 2, • • • q — 1. Then the autocorrelation 
function can be written as 

A(t) = £> A (t) (35) 

A=0 

and 

P A (t) = ^ Cfe(?+A e-^W (36) 
k 

where 

p = W fcg+A+f w 

kq+A (kq + A)\r(kq + A + l+fi)I^2yrjJl) 1 ) 

and (f>k q +A(t) is given by (34) replacing n by kq + A. For quantum numbers that are multiples 
of the revival order i.e.n = kq, we have 



<t>k q {\) = 2n(kfi + k 2 q) = 0(mod27r) 

^o(^) = Y.k C kq 

When A 7^ 0, the phases of the corresponding states are 

0fc 9+ A(J) = 27r(^A + A 2 )/q 

P A (±) = exp[-27T ? (/iAg- 1 + A 2 g- 1 )E fe c fcg+A 



(38) 



(39) 



where in obtaining (37) and (38) we have made use of the result 

2ir(k(fi + 2A) + k 2 q) = 0(mod2n) (40) 

Thus from (38) and (39) it follows that, around time t — ^ the Gazeau-Klauder coherent 
state (10) splits into q packet fractions such that the A = packet fraction has zero phase 
while relative to this the other packet fractions have a constant phase (/zAg -1 + q~ l A 2 ). 
It must be mentioned that the way the packet is divided into fractions is determined only 
by the order of the revival analyzed. Figure 4 and Figure 5 give the time dependence of 
the survival functions |Pa(^)| 2 in units t rev for A = 0, 1, 2, 3 during the whole period of the 
complete revival for fixed J = 10 and \x = 28, 80 respectively. It is to be noted that fractions 
with odd and even A differ by their revival periods. 
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Figure 4: Plot of the Survival functions \P A (t)\ 2 given in (37) for A = 0, 1,2,3 and J = 
10, li = 28. 
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Figure 5: Plot of the Survival functions \P A {t)\ 2 given in (37) for A = 0, 1,2,3 and J = 
10, n = 80. 

The intensity of the fractional revival can be defined as the value of the wave packet 
survival function at the moment of revival. This value is composed of two terms: the sum 
of survival functions for the packet fractions ( the diagonal term) and the interference term 
describing the interaction of the packet fractions 

■IH) p+ ££M;M;) 

In Figure 6(a) and Figure 6(b) we present the time dependence of the diagonal terms in 
units of t rev for fixed J = 10 and \x = 28, 80 respectively. It is seen from the figures that the 
diagonal term related to the individual packet fractions takes positive values only near the 
moments of fractional revival. 

Figure 7(a) and Figure 7(b) show the time dependence of the interference terms in units 
of t rev for fixed J = 10 and \x = 28, 80 respectively. It is seen from the figures that the 
interference term plays a constructive, destructive or indifferent role near the moments of 
fractional revival. 



i - 
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Figure 6: Plots of the diagonal term Si(t) = S^ =0 |P A (^)| 2 of the survival function for 
(a) J = 10, n = 28 and (b)J = 10, fj, = 80. 
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Figure 7: Plots of the interference term S 2 (t) = S^ =0 Ep^ A P A (^)-Pr(^) °f the survival 
function, (a) J = 10, fj, = 28 ; (b) J = 10, n = 80. 
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4 Summary 



In this paper we have constructed the coherent states for nonlinear oscillator via Gazeau- 
Klauder formalism [45,46]. These coherent states are shown to satisfy the requirements of 
continuity of labeling, resolution of unity, temporal stability and action identity. The plots 
of the weighting distribution for these coherent states are almost Gaussian in nature. The 
Mandel parameter Q is sub-Poissonian which indicates that the coherent state (10) exhibits 
squeezing for all values of /i. The fractional revivals of the coherent states are evident from 
the Figure 3(b) and 3(c) of the squared modulus of the autocorrelation function. This is in 
contrast to the result obtained in ref [84] where the wave packet revival in an infinite well for 
the Schrodinger equation with position dependent mass [84] was studied. In ref [84] it was 
found that though full revival takes place, there is no fractional revival in the usual sense. 
Instead, a very narrow wave packet is located near one wall of the well, when the mass is 
higher. In the present paper the times of appearance, the spectral compositions and the 
intensities of the fractional revivals are determined by phase analysis. 
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